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Abstract. - We present experiments on the slow growth of a single crack in a fax paper sheet 
submitted to a constant force F. We find that statistically averaged crack growth curves can 
be described by only two parameters : the mean rupture time r and a characteristic growth 
length £. We propose a model based on a thermally activated rupture process that takes into 
account the microstructure of cellulose fibers. The model is able to reproduce the shape of the 
growth curve, the dependence of £ on F as well as the effect of temperature on the rupture 
time r. We find that the length scale at which rupture occurs in this model is consistently close 
to the diameter of cellulose microfibrils. 



Introduction. - A critical stress threshold is often defined to characterize material resis- 
tance to rupture. When the applied stress is smaller than the threshold, experiments show 
that rupture can still occur after a delay in time which decreases with the applied stress and 
with temperature [1—3]. A widespread approach has been to relate time-dependent fracture 
to the creep properties of the material [4,5]. In this general framework, the creep law is an 
empirical material dependent property, the most common dependence on the applied stress 
proposed in the literature being either a power law or an Arrhenius law. More recently, several 
authors have come back to a simpler situation where the material is elastic and subcritical 
rupture is thermally activated. The focus has been mainly on prediction of the rupture time 
(or lifetime) [6-8], but has also been extended to the slow growth of a single crack in brittle 
materials [9,10]. In this Letter, we present experimental results on the slow growth of a single 
crack in a sheet of paper which behaves as a quasi-brittle material. We show that the ex- 
perimental results are compatible with a model of thermally activated crack growth in brittle 
materials when we take into account the microstructure of cellulose fibers in paper. 

Experimental set-up. Crack growth is obtained by loading in mode 1 at a constant 
force F a sheet of fax paper ( Alrey) with an initial crack in the center (fig. [1^) . The sam- 
ple dimensions are : height h = 21cm, width w = 24cm, and thickness e = 50^m. The 
experimental set-up consists of a tensile machine driven by a motor (Micro Controle U E42) 
controlled electronically to move step by step (Micro Controle ITL09). The paper sheets are 
mounted on the tensile machine with both ends attached with glue tape and rolled twice over 
rigid bars clamped on jaws. The motor controls the displacement of one jaw (400 steps per 
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micrometer) while the other jaw is rigidly fixed to a force gage (Hydrotonics-TC). The tensile 
machine is placed in a box with a humidity level stabilized at 5%. In order to work on samples 
with the same initial crack shape and length Lj, we use calibrated razor blades mounted on 
a micrometric screw and we initiate a macroscopic crack precisely at the center of the sheet. 
The samples are loaded by increasing the distance between the jaws such that the resulting 
force F is perpendicular to the initial crack direction. 




Fig. 1 - a) Sample geometry, b) Linear dependence between applied force and elongation until rupture. 

A feedback loop allows us to adjust the displacement in order to keep the applied force F 
constant with a precision better than 0.5N and a response time less than 10ms. From the area 
A of a cross-section of the sheet, A being approximatively constant, we calculate the applied 
stress a e = F/A. 

Physical properties of paper. - Sheets of fax paper break in a brittle manner. This is 
evidenced by the elastic stress-strain dependence which is quasi- linear until rupture (fig. [TJd) . 
Another sign that rupture is essentially brittle is given by the very good match between the 
two opposite lips of the fracture surfaces observed on post-mortem samples. 

A sheet of paper is a complex network of cellulose fibers. Scanning electron microscopy on 
our samples shows fiber diameters between 4 and 50/xm with an average of 18/im. Cellulose 
fibers are themselves a bundle of many microfibrils. Cellulose microfibrils have a cristallinc 
structure (therefore, they are very brittle) and are consistently found to have a diameter 
d = 2.5nm [11]. 

The mechanical properties of paper depend crucially on the humidity rate. To get re- 
producible results, the fax paper samples are kept at least one day at a low humidity level 
(< 10%) and during the experiment (~ 5%). At constant humidity level (h ~ 5%) and room 
temperature, the Young modulus of the fax paper sheets is typically Y = 3.3 10 9 N.m~ 2 . 

Direct observation and image analysis. - We light the samples from the back. A high 
resolution and high speed digital camera (Photron Ultima 1024) collects the transmitted light 
and allows us to follow the crack growth. We observe that the global deformation of the paper 
sheet during a creep experiment is correlated in a rather reproducible way to the crack growth 
whatever the rupture time. We use this property to trigger the camera at fixed increment of 
elongation (one micron) rather than at fixed increment in time. This avoids saturation of the 
onboard memory card when the crack growth is slow and makes the acquisition rate faster 
when the crack grows faster and starts to have an effect on global deformation. We acquire 2 
frames at 250fps at each trigger and obtain around one thousand images per experiment. 
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Image analysis is performed to extract the length of the crack projected on the main 
direction of propagation, i.e. perpendicular to the direction of the applied load (fig. 
Although the crack actually follows a sinuous trajectory, its projected length L gives the main 
contribution to the stress intensity factor K which we compute as: K = g(irL/2h)a e yJirL/2, 
where g(x) = ytanaT/a; takes into account the finite height h of the paper sheet [12]. 
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Fig. 2 - Extraction of the projected crack length L from the crack contour detected. 



Single crack growth. For a given initial crack length subcritical crack growth is 
obtained by applying a constant force F so that K{Li) is smaller than the critical rupture 
threshold of the material K c above which fast crack propagation would occur. During an 
experiment, the crack length increases, and so does the stress intensity factor K(L). This will 
cause the crack to accelerate until it reaches the critical length L c for which K(L C ) = K c . 

On fig. [3^ we show a typical growth curve obtained during a creep experiment with an 
applied force F = 270 N and an initial crack length Li = 1cm. Since time to rupture r is 
a statistical quantity, we prefer to plot time evolution as a function of the crack length. We 
observe that the crack growth is actually intermittent. Essentially, there are periods of rest 
during which the crack tip is pinned and does not move, and also moments when the crack 
suddenly opens and advances of a certain step size s. The crack advances by jumps until it 
reaches a critical length L c where the paper sheet breaks suddenly. Measurements of L c are 
used to estimate the critical stress intensity factor K c — 6 ± 0.5 MPa.m 1 / 2 . Beyond L c , the 
crack runs across the whole sample (about 18cm in this case) in less than one second, with 
a crack speed v > 5m.s _1 . For the same experimental conditions (same stress, same initial 
crack length, same temperature and same humidity rate), we observe a strong dispersion in 
growth curves and in lifetime while the critical length seems to be rather well defined (see 
insert in fig. [3|). In order to characterize both the average crack growth and the stepwise 
growth dynamics, a statistical analysis is required. In this Letter, we are going to focus our 
study on the average dynamics. The intermittent dynamics and in particular the step size 
statistics has already been described elsewhere [13]. 

Statistically averaged crack growth. We have performed an extensive study of crack 
growth, varying the initial crack length from Li — 1cm to Li = 4cm and the applied force 
between F — 140N and F = 280N (corresponding to an initial stress intensity factor between 
Ki = 2.7MPa.m 1 / 2 and K% — 4.2MPa.m 1 / 2 ) and repeating 5 to 20 experiments in the same 
conditions (stress, initial crack length, temperature and humidity rate). The resulting mea- 
sured lifetime varied from a few seconds to a few days depending on the value of the applied 
stress or the temperature. In order to characterize the average growth dynamics, we examine 
for given experimental conditions the average time (i)(L) the crack takes to reach a length L. 

Even though the lifetime distribution is large and the growth dynamics intermittent, the 
average growth offers a regular behavior. A typical mean growth, obtained by averaging ten 
experiments in the same conditions, is plotted in fig. [3Jd. This growth is qualitatively very 
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Fig. 3 - a) Typical stepwise growth curve for a creep experiment with an initial crack length Li = 1cm 
submitted to a constant load F — 270N. The lifetime of the sample is r = 500s and the critical length 
L c = 3.3cm. In insert, a strong dispersion is observed in crack growth profile and in lifetime for 3 
creep experiments realized in the same conditions, b) Statistical average of growth curves for 10 creep 
experiments realized in the same conditions (Li = 1cm, F — 270N). The dashed line corresponds to a 
fit using equation eq. |T} with a single free parameter £ = 0.41cm. Insert: rescaled average time (t)/r 
as function of rescaled crack length (L — £j)/C f° r various initial crack lengths and applied stresses. 
The solid line corresponds to eq. JTJ. 



close to an exponential evolution: 

(1) 

Indeed, we obtain a very good fit of the data in fig. [3Jd with eq. ([1]) setting the mean 
lifetime to the experimentally measured value and using £ as a unique free parameter. Using 
the same procedure, we extract the characteristic growth length £ for various experimental 
conditions. In the insert of fig.[3jD, rescaling the crack length by £ and the time by t for many 
different experimental conditions, we show that the data collapse on the functional form given 
by eq. (fTJ) . Moreover, we have checked that the deviation from the predicted average behaviour 
is smaller when increasing the number of experiments. 

We find that the experimental value of £ is approximatively proportional to 1 /F 2 (fig. ■ 
We also observe that the rupture time for all the experiments performed at room temperature 
is essentially a function of the initial stress intensity factor i"Q (symbols with an error bar on 
fig-Et 1 - However, for a fixed value of applied force F and initial length Li, varying temperature 
between 20°C and 120°C leads to variations of the rupture time up to four order of magnitude 
(symbols without error bar on fig. EK) • 

Model of thermally activated crack growth. - In this part, we are going to describe a model 
that is able to reproduce several of the experimental features that we have just described. 

First, we recall that the experimental crack growth curves have been obtained by perform- 
ing a statistical average on many experiments in the same condition. The model will also be a 
statistical one. We are going to consider that in average the crack growth proceeds by breaking 
individual cellulose fibers at the crack tip. Still in average, the fibers who resist to the applied 
load will be considered parallel to each other and perpendicular to the crack direction. If A is 
the typical diameter of a fiber and tl is the mean rupture time for the cellulose fiber at the 
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Fig. 4 - a) Experimental value of £ extracted from the average growth profile as a function of 1/F 2 . 
b) £ as a function of the prediction of the thermally activated rupture model. The line represents 
the best linear fit y = x/V . Its slope permits us to obtain a characteristic length scale for rupture: 
V 1/3 ~1.5nm. 



crack tip, the mean crack velocity is : v — X/tl- In that way, the increase in crack velocity 
as the crack length increases will correspond to a decrease of the rupture time tl ■ 

Whatever is the rupture mechanism, we expect that the rupture time T£ is going to depend 
on the applied load / on the fiber at the crack tip. Equivalently, we can refer to the applied 
stress a L — f /S where S ~ A 2 is the cross-sectional area of the fiber. We estimate <t l assuming 
that the fiber acts as a cohesive zone of length A and that the rest of the material behaves as 
an elastic continuum for a crack with a length L + X. Following [14], we consider that the crack 
is at equilibrium if the elastic field presents no divergence of stress. This condition imposes 
that the stress level in the fiber is controlled by the size of the fiber and can be estimated as: 
VLT\/y/\ ~ a e y/T/X. In our calculations, we will take ol — K/^/2tt\ where K is 
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Fig. 5 - a) Logarithm of the mean rupture time r as a function of K c — Ki. Data points without error 
bars correspond to non-averaged measurements obtained by varying temperature with Li = 2cm and 
various fixed values of F. b) Logarithm of lifetimes as a function of the dimensional factor AU/VkBT 
predicted by eq. ([2]| for different values of Li, F and T. The best fit log r oc AU jksT (dashed line) 
slope gives an estimation of the characteristic length scale V 1 ^ 3 ~ 2.2nm. 
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the stress intensity factor. 

A common approach to describe the creep rupture of materials is to introduce a time- 
dependent creep compliance or a rate of rupture which is a power law of the applied stress. 
Instead of assuming a phenomenological law for the creep behavior, we are going to use a 
statistical approach that has been developed in the case of brittle materials to describe the 
rupture of fiber bundles. This approach is appropriate in the case of a cellulose fiber because it 
is actually a bundle of many identical brittle microfibrils. The idea of the rupture mechanism is 
that at equilibrium there are always statistical fluctuations of stress with a variance depending 
on the actual temperature of the material. Taking into account the probability that such stress 
fluctuations exceed the rupture threshold f c of individual fibers in the bundle, it is possible 
to make some prediction on the rupture dynamics of a fiber bundle. In an homogenous fiber 
bundle where all the fibers have the same rupture thresholds, it was shown that the rupture 
time is [15, 16]: 

■(/c-/ ) 2 " 



exp 



2ak B T 



where k B is Boltzmann constant, T is temperature, a is the elastic stiffness of the fibers, /o 
is the applied force per fiber when the bundle is still intact, and f c is the rupture threshold of 
the fibers. We can rewrite this expression in term of the Young modulus Y, the applied stress 
on the fiber <jl and the rupture threshold a c . If s = nd 2 /4 is the section area of a microfibril, 
we have /o = cr^s and f c — a c s. Rupture in a microfibril is expected to occur in a small piece 
of the fiber with a length £ much smaller than the total fiber length. The relation between the 
stiffness a of this piece and the Young modulus is then: a = Ys/l. Introducing the volume 
V = si, we see that the rupture time can be written as : 



t l ~ exp 



(<x e - a L ) 2 V 
2Yk B T 



It must be understood that in this formula V represents the volume in which rupture of the 
microfibril actually occurs. Since the microfibril is broken whenever we cut a section through, 
it is physically reasonable to consider that the length t of the piece involved in rupture will 
be comparable to the microfibril diameter d. Thus, we will expect to have V ~ d 3 . 

The expression obtained for tl, gives an exponential dependence of the crack growth veloc- 
ity. If we integrate this velocity following the procedure described in [10], we recover eq. ([1]) 
where the rupture time is : 



T 



expf^, with AU=^~^ V (2) 
P \k B T ) ' 2Y y J 



and the characteristic growth length £: 

_ 2Yk B T U 

V o l (p c -o l ) 

where <ii is the value of <jl when L — Li. 

Writing explicitly the dependence of stresses on F, we obtain from eq. (J3J) that £ oc 1 /F 2 in 
agreement with the experimental observations (fig.H^). We have also £ oc yjLi/(L c — Li) but 
this dependence is difficult to verify experimentally because the ratio L c /Li does not change 
much. In fig. [4)d, we see that there is a reasonable agreement between the experimental value 
of C and the theoretical one (V — 2YkBTLi/[(Ji(a c — <7j)]. Using V as a free parameter, we 
find a characteristic scale V 1 ^ 3 — 1.5nm close to the microfibril diameter d. 
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The form of the energy barrier AU in the Arrhenius law (eq. ©) is analogous to the one 
obtained for example in [9] . Note also that it is not simply a function of the external load a e as 
in [6-8]. In fig. [5)3, we plot the rupture time as a function of AU /VksT — (cr e — <7i) 2 /ZYksT . 
We observe that the data points obtained by varying temperature are now collapsing with the 
one obtained at fixed temperature (the dispersion in this figure remains large because data 
as a function of temperature and constant F and Li have been averaged less than those at 
constant temperature for technical reasons). From a fit of the data, we obtain independently a 
new estimate of V which gives a characteristic scale V 1 ^ ~ 2.2nm. Once again this estimate 
is close to the microfibril diameter d. 

Conclusion. - We have studied the subcritical growth of a single crack in a sheet of paper 
during creep experiments. We have shown that statistically the crack growth curves can be 
well described with only two parameters: the rupture time r and a characteristic growth 
length £. We have proposed a model which is able to reproduce the shape of the growth 
curves as well as the dependence of C with the applied load and the influence of temperature 
on the rupture time. The model takes advantage of the fact cellulose fibers are actually a 
bundle of many small brittle fibers with a diameter of nanometer scale. A previous prediction 
of the rupture time for thermally activated rupture of brittle fibers has been used to derive the 
crack velocity. We find that the experimental data are compatible with the idea that rupture 
occurs at a nanometric scale in the cellulose fibers. While the model we used work rather 
well, it remains to understand if the experimental results could be as well described by time- 
dependent fracture models based on the material rheological properties. Such a study would 
require for example a detailed investigation of the physical properties of individual fibers. 

* * * 
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